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Abstract 

We argue that the Einstein-Yang-Mills theory presents nontrivial solutions with a NUT charge. These 
solutions approach asymptotically the Taub-NUT spacetime. They are characterized by the NUT param¬ 
eter, the mass and the node numbers of the magnetic potential and present both electric and magnetic 
potentials. 

The existence of nontrivial Einstein-Yang-Mills solutions with NUT charge in the presence of a nega¬ 
tive cosmological constant is also discussed. We use the counterterm subtraction method to calculate the 
boundary energy-momentum tensor and the mass of these configurations. Also, dyon black hole solutions 
with nonspherical event horizon topology are shown to exist for a negative cosmological constant. 


1 Introduction 

According to the so called ” no-hair” conjecture, an asymptotically flat, stationary black hole is uniquely 
described in terms of a small set of asymptotically measurable quantities. It was believed that this property 
of black holes, following from the uniqueness and nonexistence results proven rigorously for a variety of 
field theories, persists for general matter sources. The hypothesis was disproved more than ten years ago, 
when several authors presented a counterexample within the framework of SU(2) Einstein-Yang-Mills (EYM) 
theory 0. Although the new solution was static and had vanishing Yang-Mills (YM) charges, it was different 
from the Schwarzschild black hole and, therefore, not characterized by its total mass. Despite its lack of 
stability, this nonabelian solution still presents a challenge to the no hair conjecture. Since then, a lot of 
black hole with different hairs have been found (for a general review of such solutions see 00). 

It is however worth inquiring, what happens if we drop asymptotic flatness? Will these ” hairy” black 
holes survive? As proven by various authors, this is the case in the presence of a cosmological constant A. In 
the last years, nonabelian black hole solutions in a cosmological context were considered by various authors. 
Remarkably enough, for an asymptotically anti-de Sitter (AdS) spacetime, there are even stable solutions 

BI¬ 
AS discussed in 0 0 , a variety of well known features of asymptotically self-gravitating nonabelian 
solutions are not shared by their AdS counterparts. For example, there are finite energy solutions for several 
intervals of the shooting parameter (the value of gauge function at the event horizon), rather than discrete 
values. There are also stable monopole solutions in which the gauge field has no zeros, and the asymptotic 
value of the gauge potential is not fixed. This behavior is drastically different from those observed for EYM 
black holes in asymptotically flat I] or de Sitter spacetime 0. Their regular counterparts are discussed in 
00 and present also interesting properties. 

In some sense, the minimal deviation from the asymptotically flat case is to include a ’’dual” or ’’magnetic” 
mass in the theory. That general relativity does admit configurations with ’’magnetic” masses has been clear 
after the discovery of the NUT solution 00 (see also nni for a recent discussion of this solution). This metric 
has become renowned for being ” a counterexample to almost anything” El and represents a generalization 
of the Schwarzschild solution El It is usually interpreted as describing a gravitational dyon with both 
ordinary and magnetic mass. It is the NUT charge which plays a dual role to ordinary mass, in the same 
way that electric and magnetic charges are dual within Maxwell theory El As discussed by many authors, 
the presence of magnetic-type mass (the NUT parameter n) introduces a ’’Dirac-string singularity” in the 
metric (but no curvature singularity). This can be removed by appropriate identifications and changes in 
the topology of the spacetime manifold, which imply a periodic time coordinate. The Killing symmetries of 
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this solution are still translations and SO(3) rotations. Spherical symmetry in a conventional sense is lost, 
though, when the NUT parameter is nonzero, since the rotations act on the time coordinate as well. The 
periodicity of the time coordinate prevents an interpretation of NUT metric as a black hole. Moreover, the 
metric is not asymptotically flat in the usual sense although it does obey the required fall-off conditions. 

An alternative suggestion due to Bonnor is to regard the string as a physical source HI- In this approach, 
the source of the field is a spherically symmetric mass together with a semi-infinite massless source of angular 
momentum along the symmetry axis. 

A large number of papers have been written investigating the properties of the gravitational analogs of 
magnetic monopoles. The NUT solution has been generalized in a number of ways m A charged solution 
of the Einstein-Maxwell equations with magnetic mass has been found by Brill m • There are also nontrivial 
solutions of the low-energy string theory with NUT charge (see e.g. (I7J). I 11 HI a general framework was 
formulated to handle spacetimes with dual mass, and showed that the dual mass is absolutely conserved in 
the presence of gravitational radiation. Unfortunately, the pathology of closed timelike curves (CTC) is not 
special to the NUT solution but afflicts gravitation solutions with ’’dual” mass in general m The existence 
of magnetic mass requires either that the metric (not the curvature) be singular on a two-dimensional timelike 
world sheet or that the spacetime contains CTC. 

As emphasized in m, this condition emerges only from the asymptotic form of the fields, and is com¬ 
pletely insensitive to the precise details of the nature of the source, or the precise nature of the theory of 
gravity at short distances where general relativity may be expected to break down. The same construction 
can be made for a spacetime with a negative cosmological constant A, since it has also a well defined spacelike 
infinity generated by a timelike Killing vector. A similar analysis will yield the same pathologies. 

Consequently, spacetimes with NUT charge are not likely to be of interest as models for macroscopic 
objects. Nevertheless, there are various features suggesting that their role cannot be neglected in the context 
of quantum gravity m 

A natural question arises: do there exist counterparts of the well-known hairy black hole solutions 
possessing a magnetic mass? And what new effects emerge due to the presence of a NUT charge? 

Solutions of the YM equations on the Euclidean section of the Taub-NUT spacetime have been considered 
in the literature Ea, while a study of their Lorentzian counterparts is still missing. This fact is due perhaps 
to the widely belief that these spacetime are unphysical and thus cannot accomplish for interesting effects. 

However, at the moment there are many indications that the NUT charge is an important ingredient 
in low energy string theory El, conclusion enhanced by the discovery of ” duality” transformations which 
relate superficially very different configurations. In many situations, if the NUT charge is not included in 
the study, some symmetries of the system remain unnoticed (see e.g. m)- Therefore, we may expect to find 
a similar picture when discussing duality properties (yet to be discovered) of gravitating nonabelian fields. 

There are also a number of reason to extend this analysis to configurations with a negative cosmological 
constant. Black holes with hair in such theories are useful for probing not only quantum gravity, but also 
may be interesting tests of AdS/CFT, particularly for classical stable configurations. Also, when A < 0, 
the so-called topological black holes, whose topology of the event horizon is no longer the two-sphere S 2 
may appear (see m for reviews of the subject). These black holes embedded in ’locally AdS’ background 
spacetimes (background locally isometric to spacetimes of constant negative curvature) have been seminal to 
recent developments in black hole physics. Vacuum Taub-NUT-AdS solution in four dimensions provided the 
first test bed for AdS/CFT correspondence in spacetimes where the asymptotic structure was only locally 
asymptotic AdS mm- Their Euclideanized form have been a subject of interest in recent years, in part 
because these asymptotically locally AdS spaces are relevant to 2+1 dimensional ’’exotic” CFT that live on 
the world volume of M2 branes after placing them on a squashed three sphere EH!. It is therefore natural 
to suppose that generalization of these spacetimes for a nonabelian matter content might provide us with a 
better understanding of these issues. 

In this paper we restrict attention to SU( 2) EYM theory with A < 0, and present numerical arguments for 
the existence of nontrivial solutions with NUT charge. For a negative cosmological constant, these solutions 
present nonvanishing nonabelian charges (with a suitable definition of charges in this case). We find a family 
of solutions generalizing for a NUT charge the known spherically symmetric dyon black hole solutions 0. 
The Taub-NUT-AdS family of metrics contains also solutions where the angular spheres (9, ip) are replaced 
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by planes or hyperboloids. For vanishing NUT charge, these solutions correspond to topological black holes. 
In this paper we construct numerically the nonabelian dyon counterparts of these configurations, generalizing 
for a NUT charge the topologically nonabelian black holes found in m ■ The properties of these solutions do 
not differ too much from the properties of the dyon black hole solutions (with n = 0) in AdS spacetime. In 
this case we find again a continuum of solutions possessing nonvanishing nonabelian charges and arbitrary 
mass. 

However, the configurations with A = 0 and NUT charge preserve the basic properties of the well known 
hairy black hole solutions. The solutions do not posses nonabelian charges (although a nonzero electric 
potential is always present) and, for a given n, the solution is uniquely determined by the mass and the node 
number of the magnetic potential. 

The paper is structured as follows: in the next Section we explain the model and derive the basic 
equations and boundary conditions. In Section 3 we show results obtained by numerical calculations and 
discuss the general properties of solutions as well as particular cases. The boundary stress tensor and the 
associated conserved charge are computed in Section 4, by using a counterterm prescription. We discuss also 
the asymptotics of these spacetimes and find that a (2 + l)-dimensional Godel-type universe is the boundary 
of a topologically NUT-AdS spacetime. We conclude with Section 5 where the results are compiled. 


2 General framework and equations of motion 

2.1 Metric ansatz and symmetries 

The EYM coupled system is described by the action 


J = / d A Xy/^g 


—-—(R - 2A) - -F°;„F atlv 
16ttG v ’ 4 ^ 


f d 3 x\/~-hK. 
87rG Jq M 


(i) 


Here G is the gravitational constant, R is the Ricci scalar associated with the spacetime metric g and 
F IJV = \r a F- v is the gauge field strength tensor defined as 

F/J.V = d^Av — d v A^ + ie [A^, A u ] , (2) 

where the gauge field is A M = fr a A“, r a being the Pauli matrices. The last term in 0 is the Hawking- 
Gibbons surface term 129 . where K is the trace of the extrinsic curvature for the boundary dAd and h is 
the induced metric of the boundary. Of course, this term does not affect the equations of motion but is 
important when discussing quantum properties of gravitating solutions. 

Variation of the action 0 with respect to the metric g^ leads to the Einstein equations 

R^ - A g v = 8 ttGT^, (3) 

where the YM stress-energy tensor is 

TV = 2T \{F^F v gg^ - -^g^F^F^). (4) 

Variation with respect to the gauge field A^ leads to the matter field equations 


+ ie[A^ EH = 0. 


(5) 


We consider spacetimes whose metric can be written locally in the form 

dr 2 f) 

ds2 = N^) + R2 ^ d ° 2 + A 2 ( 0 )^ 2 ) - N(r)e~ 25 ^(dt + 4 nf 2 k {-)d^) 2 . (6) 
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The discrete parameter k takes the values 1, 0 and —1 and implies the form of the function fk{9) 


fm 


sin 0, for k = 1 
9 , for k = 0 

sinh0, for k = — 1. 


( 7 ) 


We define the NUT parameter n, in terms of the coefficient appearing in the differential dt + Anf 2 (9/2)dtp. 

For n = 0, metrics of this form have been considered by various authors when discussing topological 
black holes with A < 0. In this case, dH'jU = d0 2 + f 2 {9)d<p 2 is the metric on a two-dimensional surface E 
of constant curvature 2k. Thus, without a NUT charge, the topology of spacetime is R 2 x H 2 , where H 2 is 
the topology of the surface E. When k = +1, the universe takes on the familiar spherically symmetric form, 
and the ( 9 , <p) sector has constant positive curvature. When k = 0, the E is a flat surface, while for k = — 1 
the ( 9 , ip) sector is a space with constant negative curvature, also known as a hyperbolic plane. A discussion 
of these cases is available in m, When E is closed (which we’ll suppose to be always the case), we denote 
its area by V, defining V = f fk(9)d9dp. 

The form of N(r), S(r) and R(r) emerges as result of demanding the metric to be a solution of the EYM 
equations with negative cosmological constant (note the existence of a metric gauge freedom which can be 
used to fix the expression of R{r)). 

For this ansatz, the nonvanishing components of the Einstein tensor E” read 
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where a prime denote a derivative with respect to r. 

Apart from the Killing vector dt , this line element possesses three more Killing vectors characterizing the 
NUT symmetries 


Ki 

■ a , /*(<?)« 10 fk(B/2) a 

= sin yd e + cos ip / d<p + 2 n cos <p d t , 

Jk\“) 


I<2 

a ■ f' k {0 ) n 0 . f k (9/2) a 

= cos (pdo sin ip dp 2nsm<p d u 

Jk{“) 

(9) 

K-i 

= dp — 2nd t . 



These Killing vectors form a subgroup with the same structure constants that are obeyed by spherically 
symmetric solutions [K a ,Kb] = e a b c K c - 

The properties of case k = 1 are well known from the study of Taub-NUT solution. The rotational 
symmetries act in the usual way on the angular coordinates, but also involve time translations in order 
to preserve the differential dt + An sin 2 (0/2)d(^. Thus we have lost spherical symmetry, at least in the 
conventional sense. 

Again, the n sin 2 (6/2) term in the metric means that a small loop around the z—axis does not shrink 
to zero at 9 = 0 and at 9 = ir. These singularities can be regarded ar the analogue of a Dirac string in 
electrodynamics and are not the usual degeneracies of spherical coordinates on the two-sphere. This problem 
was first encountered in the vacuum NUT metric. One way to deal with this singularity has been proposed 
by Misner m- His argument holds also independently of the precise functional form of N and <5. In his 
construction, one considers one coordinate patch in which the string runs off to infinity along the north axis. 
A new coordinate system can then be found with the string running off to infinity along the south axis. This 
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indicates that the string is an artifact resulting from a poor choice of coordinates. In the intersection of 
these two coordinate patches, a necessary condition for consistency is to identify the time coordinate with 
period 8nn. The coordinate singularity at 6 = it can be removed via the introduction of 

t' = t + 4rup. ( 10 ) 

It is clear that the t coordinate is also periodic with period 87 m and essentially becomes an Euler angle 
coordinate on S 3 . Thus there are CTC which can not be removed by going to a covering space and no 
reasonable spacelike surface. One finds also that surfaces of constant radius have the topology of a three- 
sphere, in which there is a Hopf fibration of the S 1 of time over the spatial s 2 m- 

Therefore for n different from zero, the k = 1 metric structure © shares the same troubles exhibited 
by the Taub-NUT gravitational field in the presence, if any, of an electromagnetic field, i.e. the solutions 
cannot be interpreted properly as black hole. However, for k = 0,-1, the fibration is trivial and there are 
no Misner strings m, but again, there are no reasonable spacelikc surfaces. 

One can remark another interesting property of the line element ©• The static nature of a hairy 
spherically symmetric black hole solution implies that it can only produce a ” gravitoelectric” field. In 
the case of a NUT charge, the existence of the cross metric term dpdt shows that this space has also a 
’’gravomagnetic” field. Thus the line-element © it is Kerr-like in the sense that it has a crossed metric 
component g v t which generates ’’gravomagnetic” effects. This term does not produce an ergoregion but it 
does produce an effect similar to the dragging of inertial frames 1321 . 


2.2 577(2) connection 

The most general expression for the appropriate SU( 2) connection is obtained by using the standard rule for 
calculating the gauge potentials for any spacetime group mm- According to Forgacs and Manton, a gauge 
field admit a spacetime symmetry if the spacetime transformation of the potential can be compensated by 
a gauge transformation ESi. 

Taking into account the symmetries of the line element © we find the general form 


A = 


2e 


+ 4nfk(^)dtp)u(r)T 3 + v{r)r 3 dr + (w(r)n + ui(r)T 2 ^d6 + fk(0)r 3 + (w(r)r 2 - w(r)ri)/ fe ( 6 >) dtp j , ( 11 ) 


(see Ref. jS| for an extensive discussion of this ansatz in the case k = 1, n = 0). The gauge connection 
CD remains invariant under a residual f/(l) gauge symmetry which can be used to set v = 0. Because 
the variables to and to appear completely symmetrically in the EYM system, the two amplitudes must be 
proportional and we can always set Co = 0 (after a suitable gauge transformation). By using the residual 
U{ 1) gauge symmetry, we can also consistently set v = 0. Hence, without any loss of generality, the gauge 
potential is described by two functions to(r) and u(r ) which we shall refer to as magnetic and electric 
potential, respectively. 

Therefore, the components of the field strength tensor are 

Fro = ^ = yw'A(fl) + |W/ k 2 (^ F rt = ^u', (12) 

Fs^p = -y4 nutofki^) + y(w 2 - k + 2 nu)f k (0), F et = -yww, F vt = yuw/ fc (0). 

It is difficult to find a satisfactory definition of the nonabelian electric and magnetic charges in the presence 
of a NUT parameter. For n = 0, the usual (-gauge dependent) definition is 



where F^ v is the dual field strength tensor. However, the integral G3J does not intrinsically characterize 
a field configuration, as it is gauge dependent 0. The definition of conserved currents and charges in a 
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non-Abelian Yang-Mills theory is a problem approached by different authors in the last decades (see e.g. 
G2H2E]), with various solutions. A gauge invariant definition for the nonabelian charges has been proposed 
in [37] (see also m) 



d9dtp\F S(p \, Q m 


e 

47r 


d9dip\F eip \, 


(14) 


where the vertical bars denote the Lie-algebra norm and the integrals are evaluated as r —■> oo. One can verify 
that this definition yields the absolut value of the charges as defined in m- Also, for the field configuration 
we consider here, the same expression G3 for Qe, Qm is obtained when making use of the existence of a 
Lie-algebra valued ” Killing scalar” [2] SHI ■ 

In the presence of a dual mass (n ^ 0) it seems that we can only define the magnetic and electric charges 
in terms of the asymptotic behavior of the gauge field by analogy to that in asymptotically flat spacetime, 
i.e. from ltl2ll ~ Qe/v 2 and ~ QmA-( 0) (a similar problem occurs for an 17(1) field fTj). 

We now remove the dependence on the coupling constants G and e from the differential equations by 
using the following rescaling r —> ry/4nG/e, n —* ny/4rG/e, u —> ue/y/AirG , A —> Ae 2 /(47rG). 

For this ansatz, the nonvanishing components of the energy-momentum tensor read 
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(15) 


2.3 Field equations 

One explicit solution of the EYM equations is well known. This is the charged generalization of the NUT 
solution found by Brill m, generalized for a nonvanishing A m- These abelian configurations have a mass 
M, electric charge Q and magnetic charge P = k 


u = 0, 


nP-Qr 

u = —^, o = 0, 


r z + n 


2 ’ 


. rr. -y .. , 2 (Mr + kn 2 ) + (Q 2 + P 2 ) A r 4 + 6r 2 n 2 - 3n 4 

R{r) = V r 2 +n~, N = k - 7 U —2 -T- 2 U —2 - 

r z + n z 3 r z + n z 


(16) 


Some features of these Einstein-Maxwell solutions are discussed in ESI- The properties of the vacuum 
topologically-Taub-NUT-AdS solutions (on the Euclidean section) are discussed in ESI- 

The Brill solution is written in a metric gauge <5 = 0. For a nonabelian configuration, a more convenient 
gauge, used also in the numerical computation is the Schwarzschild one R(r) = r. For this choice, we find 
from © that the nonabelian potentials should satisfy the equations 


{Ne~ s oj')' = 
(■ e s u'r 2 )' = 

while the Einstein equations © imply 
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Note that, for n / 0, these equations do not present the usual symmetry u —► —it, which is related to the 
absence of time reflection symmetry in the metric ©■ Also, for an asymptotically flat hairy black hole, 
the vanishing of the electric potential is imposed by finite energy requirements ED This conclusion remain 
valid in the presence of a positive cosmological constant 0- However, for n / 0 we cannot consistently set 
u = 0 for nonabelian solutions. Thus, all solutions with NUT charge necessarily present an electric part and 
possibly an electric charge. 

2.4 Boundary conditions 

We want the metric © to describe a nonsingular, asymptotically NUT spacetime outside an horizon located 
at r = rv Here N(rh) = 0 is only a coordinate singularity where all curvature invariants are finite. A 
nonsingular extension across this null surface can be found just as at the event horizon of a black hole. If the 
time is chosen to be periodic, as discussed above, this surface would not be a global event horizon, although 
it would still be an apparent horizon. The parameter n describes the departure of the metric © from the 
static black hole solutions. The regularity assumption implies that all curvature invariants at r = ru are 
finite. Also, the energy density T t * is assumed to be regular everywhere. 

We find that the general properties and asymptotic form of the solutions crucially depend on the presence 
of A. 


2.4.1 The case A < 0 


For A < 0, we find the following asymptotic form of the solution 
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( 20 ) 


and ujq, uo , (Iq are constants to be determined. Thus, the geometry approaches asymptotically the 
(topological-) Taub-NUT-AdS spacetime. 

For these asymptotics we define the nonabelian charges 


(® E \ = }L( ; Si ^. (2i) 

\ Qm ) 47t \ k - Wq - 2 nu 0 J ' K 

In the rest of this paper Qe, Qm a re expressed in units of V/4tt. 

The corresponding expansion as r —> is 

N(r) = N h (r - m) + O ((r - r^) 2 ) , 
e -<5(U _ e -d 0 _|_ _ Vh j _|_ q (( r _ r/i ) 2 ) ; 

w(r) = u> h + wi(r - r h ) + O ((r - r h ) 2 ) , (22) 

u{r) = rti(r - r h ) + u 2 (r - r h ) 2 + O ((r - r h ) 3 ) , 


7 




where 


N h 

di 


Wi 

U2 


l (k _ Ar I - fcLS! _ , 


rh 

2 

— ( 


-do 


2__^1 2 d 0 ^ 

2rJ + ' Nl 


Uh(u 2 h - k) 

rlNh ’ 


“h U 1 „-2dn n K- fc ) , „-dn U l d l 

t- I it- ^ 


^r 2 


( 23 ) 


In the numerical computation we have set do = 1, without loss of generality. The usual condition <5(oo) = 0 
is connected with our choice by a a simple time-coordinate rescaling, together with a rescaling of u and n 
(we observe that the field equations m - m remain unchanged if we take t —> t/c, d —> <5 — logc, u —> 
cu, n —> n/c). The magnetic charge is not affected by this rescaling, while for the electric charge we have 
Qe —*■ cQe- 


The function N should also satisfy the obvious condition 
bound on u>h 


dN(r) 

dr 
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> 0. For given (r>j, A), this places a 
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(24) 


and implies the positiveness of the quantity Lo'{rh). In the k = — 1 case, this implies the existence of a 
minimal value of |A|, i.e. for a given rh 

|A|>-l(l + -t). (25) 

' h 1 h 

The asymptotic behavior m suggest to adopt the following form for the metric function N(r) 


N( 7’) = k — n 2 A — 
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For n = 0, m(r) corresponds to the local mass-energy function and its limit as r —> cx) corresponds to ADM 
mass. We present arguments in Section 4 that this interpretation is still valid for a nonvanishing NUT 
charge. The numerical results have been obtained for this form of N(r) which factorizes the divergent part 
Ar 2 /3. 


2.4.2 The case A = 0 


As expected, for A = 0, the behavior of the solutions at infinity is very different as compared to the case 
A < 0. The only allowed configuration for a vanishing cosmological constant is k = 1. In this case, we find 
the asymptotic form of the solutions 
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where ui \, 112 , g ?4 are constants to be determined. Again M is to be interpreted as the total mass of the 
solutions. The boundary conditions on the event horizon are obtained from (12211 by taking the limit A = 0. 

In asymptotically Minkowski spacetime (n = 0), the electric components of the YM connection are 
forbidden in static solutions El We find that this is not valid for solutions possessing a NUT parameter. 
Moreover, as seen from the YM equations, the existence of a nonzero electric potential is required in this 
case. This has to be associated with the Kerr-like form of the line-element ©■ A similar property is shared 
by axially symmetric rotating black holes in asymptotically flat spacetime ESI- 

We remark also that, for n 7 ^ 0, the decay of u is too rapid to give a net electric charge and w(oo) 2 = 1 
corresponds to zero magnetic charge. Thus, these nonabelian solutions are not dyons in the usual sense but 
do possess both electric and magnetic fields. A similar property has been noticed for the asymptotically flat 
solutions of an EYM-massive axion theory fA‘2\ . 


3 Numerical results 

Although an analytic or approximate solution still appears to be intractable, we present in this section 
numerical arguments that the known EYM black hole solutions can be extended to include a NUT parameter. 
The properties of the solutions in the cases A < 0 and A = 0 are very different and we will discuss them 
separately. 

3.1 Charged topological black holes 

We start by discussing the particular case n = 0, A < 0. The properties of the spherically symmetric dyon 
solutions k = 1 are studied in Ref. 3] (see also El)- In the presence of a negative cosmological constant, 
the ’’nonabelian baldness” theorems (forbidding SU( 2) dyons) are no longer valid: there are essentially 
nonabelian solutions with an electric sector. Dyon solutions of the EYM equations are determined if the 
parameter -ui in is chosen to be nonzero for a given A. One finds dyon solutions with finite mass for 
continuous varying u>h (and the corresponding parameter in the electric sector) and any number of nodes p 
including p = 0. 

We shall therefore focus on the cases k = 0 and k = —1. Although predicted in m , these topological 
dyon black holes are not discussed in literature. An analysis of the properties of corresponding monopole 
solutions is presented in m ■ These are black hole solutions with nonspherical event horizon topology and 
possessing a nonvanishing magnetic potential w(r). For u = 0 only nodeless solutions are allowed. Also all 
k = 0 and the k = — 1 solutions with |wo| > 1 are found to be stable against linear fluctuations. Another 
interesting property is the existence of k = —1 solutions with M < 0, a common property, however, of this 
class of topological black holes EU Black hole solutions with nonspherical topology do not have globally 
regular counterparts. 

For n = 0 and k = 0, — 1, we numerically solve the full set of Eqs. m-m with the boundary conditions 
(PI) using a standard shooting method. Starting with the horizon expansion the equations were 

integrated for a range of values of (rh, A) and varying (u>h, iti). Since in this case the field equations are 
invariant under the transformation to —> —w, u —> —u, only values of u>h, u\ greater than zero are considered. 
For sufficiently small u>h, all field variables remain close to their values for the abelian configuration with the 
same rh ■ Significant differences occur for large enough values of LOh and the effect of the nonabelian field on 
the geometry becomes more and more pronounced. 

The solutions we obtain have many common properties with their spherically symmetric A < 0 coun¬ 
terparts. Solutions are classified by the parameters M, Qm and Qe ■ The behavior of metric functions m 
and 5 is similar to the spherically symmetric monopole black hole solutions. For every considered value of 
(A, rh), we find regular black hole solutions for compact intervals of LOh and U \. For u>h > ^and u\ > u\ 
the solutions blow up or the function N(r) becomes negative. The values of u1 increase as |A| increases 

and depend also on the value of rh- There are also solutions for which ujoo > 1 although LOh < 1. 

In contrast to the monopole case, both nodeless solutions and solutions where lo crosses the axis can exist. 
This is not a surprise, since the simple argument presented in m for the existence of nodeless solutions only 
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is invalidated by the presence of an electric potential. For k = 0 we have always m(r) > 0 and the black 
holes therefore only occur with positive values of mass. Typical solutions in these cases are presented in 
Figure 1. In Figure 2 we present the mass parameter M, the nonabelian charges Qe , Qm and the values of 
the gauge potentials at infinity u o,u>q as a function of the magnetic potential value on the horizon for fixed 
T] j, A. Not unexpectedly, when the surface S has a negative curvature, we have found for suitable values of 
(A, r h , u h , m), a class of solutions with M < 0 or even M = 0, apart from configurations with positive M. 

3.2 Solutions with NUT charge and negative cosmological constant 

Our approach for a nonvanishing n is very similar to the black hole case. By using a standard shooting 
method, we integrate the field equations (TTTI-ITTTl) in the general case. The NUT charge n, as well as the 
cosmological constant A are primary input parameters. We call rh, Uh, u\ secondary input parameters. A 
complete analysis of the complex correlation between the parameters of the theory is beyond the purposes 
of this paper. Here we present only the results of a preliminary investigation. However, we found nontrivial 
solutions for a large range of primary parameters n , A, some fixed values of rh and compact intervals of 
secondary parameters uih, u±. Nontrivial solutions seems to exist for every choice of n, A satisfying P|) . 

The numerically solutions we find in this case retains all the features of the ’’gravitoelectric” case n = 0. 
Also, the fundamental properties of the solutions do not depend on the value of k. For small values of the 
NUT parameter the solutions do not considerably deviate from the corresponding black hole solutions. 

For fixed primary input parameters and a given rh, solutions are found for continuous sets of secondary 
parameters Uh and u±. Again we find a continuum of solutions where w crosses the axis an arbitrary number 
of times depending on u>h and u\. The existence of critical secondary parameters and itf is also noticed. 
Similarly there exist solutions where ui does not cross the axis. For k = — 1 we find a class of solutions with 
M < 0 or even M = 0, apart from configurations with positive M. Typical solutions are presented in Figure 
3. We can see that the electric potential starts at zero and asymptotically approaches a finite value uq. The 
mass parameter stay also finite. In Figure 4 we present the results of the numerical integration for a varying 
NUT charge, while keeping (A, ?>,, ujh, u\) fixed. In Figure 5 we vary the secondary parameter u\ for fixed 
(A, n, rh,0Jh)- Solutions with nonavanishing nonabelian charges are found also for u\ = 0 , u>h ^ 0. 

3.3 Solutions with NUT charge and A = 0 

The case A = 0, k = 1 is rather special. Such a spacetime is not asymptotically flat in the usual sense, 
although it does obey the required fall-off conditions. 

The solutions we find are a kind of extension of the well known asymptotically flat black hole solutions 
[T]. As a new feature, they present a nonvanishing electric potential but no nonabelian charges. 

We find that, given (n, rh), solutions may exist for a discrete set of shooting parameters u>q and u\. By 
using a standard ordinary differential equation solver, we evaluate the initial conditions at rh for global 
tolerance 10 —12 , adjusting for fixed shooting parameters and integrating towards r —> oo. The integration 
stops when the asymptotic limit El) is reached. 

The profile of typical solutions is presented in Figure 6. We observe that the electric potential starts at 
zero and monotonically increases rapidly to some maximal value, afterward decreasing slowly to zero. The 
behavior of the u, m, S in the region r > rh is qualitatively similar to that for the black hole solutions. 
The function u> starts from some value 0 < wj, < 1 at the horizon and after p oscillations around zero 
tends asymptotically to (—l) p . We have found solutions up to p = 3, but we have no reason to doubt the 
existence of similar solutions for every value of p. For a given n, the geometry in the asymptotic region 
r —> oo is Taub-NUT, with a mass M = M p (rh)- Similar to the black hole case, the node number p cannot 
be associated with any kind of YM charge. 

For all solutions, we found that always |w| < 1 everywhere outside the horizon. We noticed that, for a 
given p, the value of the mass-parameter M for n ^ 0, is always smaller that the corresponding black hole 
ADM mass. For example, the known value for n = 0 is M = 0.937; for n = 0.284 we found M = 0.887 while 
for n = 0.493, M = 0.754 (this is valid for p = 1, rh = 1). 
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Another point which deserves further research is the existence in this case of a maximal value of n. For 
a given value of rh we failed to find numerical solutions with the right asymptotics beyond a value n max of 
the NUT parameter. Since we deal with a two-dimensional shooting problem, it is difficult to find accurate 
values for the n max . For example, for configurations with rh = 1, p = 2, we find no regular solutions beyond 
n > 0.21. 


4 A computation of the stress tensor and the mass 

In order to compute quantities like the ADM mass one usually encounters infrared divergences, which are 
regularized by subtracting a suitable chosen background Enum Such a procedure, however, in general 
is not unique; in some cases (such as solutions with NUT charge) the choice of reference background is 
ambiguous. Recently, in order to regularize such divergences, a different procedure has been proposed mi 
This technique was inspired by AdS/CFT correspondence and consists in adding suitable counterterms I c t 
to the action. These counterterms are built up with curvature invariants of a boundary dM (which is sent 
to infinity after the integration) and thus obviously they do not alter the bulk equations of motion. The 
following counterterms are sufficient to cancel divergences in four dimensions, for vacuum solutions with a 
negative cosmological constant 


7ct 8 tt 


d 3 xV—h 


’dM 


2 l 
-l + 2 U 


(28) 


Here 1Z is the Ricci scalar for the boundary metric h. In this section we will define also l 2 = —3/A. 

Using these counterterms one can construct a divergence-free stress tensor from the total action 
/=4uik+4urf+/ct by defining 

2 SjI 1 2 

T a b = =—{K ab —Kh ab —-h ab + lE ab ), (29) 

where E ab is the Einstein tensor of the intrinsic metric h ab - The efficiency of this approach has been 
demonstrated in a broad range of examples. The counterterm subtraction method has been developed on 
its own interest and applications. If there are matter fields on M additional counterterms may be needed 
to regulate the action. However, we find that for a SU{2) nonabelian matter content the prescription (12811 
removes all divergences. The results we find by using the asymptotic expressions m for the boundary stress 
tensor at large r are 


Tbs 

Tipip 

T<pt 

T tt 


8ir r \ r 2 ) ' 

T8m„a^)I + o(-I), 



1 2M ^ 
8^17 + 0 



(30) 


Direct computation shows that this stress tensor is traceless. This result is expected from the AdS/CFT 
correspondence, since even dimensional bulk theories with A < 0 are dual to odd dimensional CFTs which 
have a vanishing trace anomaly. 

We can use this formalism to assign a mass to our EYM solutions by writing the boundary metric in an 
ADM form j2| 

h^dx»dx v = - N%dt 2 + cr ab {dx a + N“dt){dx b + N b dt). (31) 
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If V' is a Killing vector generating an isometry of the boundary geometry, there should be an associated 
conserved charge. The conserved charge associated with time translation is the mass of spacetime 


M = / dP-'-xy/oNvz. (32) 

JdY. 

Here e = n fJ "n v T fil , is the proper energy density while is a timelike unit normal to E. By using this 
relation, we can find the mass of our solutions, which is always MV/An, as expected. This result keeps valid 
in the limit k = 1, A = 0. 

The metric restricted to the boundary h ab diverges due to an infinite conformal factor r 2 /l 2 . The 
background metric upon which the dual field theory resides is 

l 2 

lab = lim -^h ab . (33) 

r— xx) r 


For the asymptotically topologically-Taub-NUT-AdS solutions considered here, the boundary metric is 

l ab dx a dx b = l 2 (dd 2 + fl(d)dy 2 ) - (dt + 4nfg(^)dtp) 2 . (34) 

In the light of AdS/CFT correspondence, Balasumbramanian and Kraus have interpreted Eq. 121111 as 
r ab = -j= where r ab is the expectation value of the CFT stress tensor. Corresponding to the boundary 
metric Gi, the stress-energy tensor r ab for the dual theory can be calculated using the following relation 

m 

\F I ll ah Tbc = lim \f~-hh ab T bc . (35) 


By using this prescription, we find that the stress tensor of the field theory has this simple form 

Tab = ^[Su a u b +i ab ], (36) 

where u a = 5/. This is the standard form for the stress tensor of a (2+1) dimensional CFT. This tensor is 
covariantly conserved and manifestly traceless. 

For the k = 1 vacuum Euclidean case, it was conjectured that the strongly coupled dual CFT is the 
infrared limit of large N, D = 3 super-Yang-Mills theory, but now defined on a squashed S 3 rather than 
the round one EElESj- However, for the EYM action m we do not know the underlying boundary CFT. In 
particular, we do not know what the SU (2) field corresponds to in CFT language. The bulk YM fields have 
nothing to do, of course, with the nonabelian fields of the dual gauge theory. 

We end this section with several considerations on the boundary line-element 03, not directly related 
to the aim of this paper. 

After the redefinitions 9 = p/l, l = A/m, n = Sl/m 2 , the line element (13411 can be written in the form 

da 2 = dp 2 + fli.mp)d<p 2 - (^fk(^-)d(p + dt) , (37) 

which is the standard form of the nontrivial (2+l)-dimensional part of an homogeneous Godel-type universe. 
This geometry can also be interpreted as resulting from the squashing of three dimensional anti-de Sitter 
(AdS) lightcones [H>| . 

The general form for a Godel-type homogeneous spacetime was found in 1983 by Rebougas and Tiomno 
m- It can be expressed as the direct Riemannian sum dz 2 + dcr 2 of a flat factor and the (2 + l)-dimensional 
metric 03- The famous Godel rotating universe @3 corresponds to the case k = — 1, to 2 = 2H 2 (l 2 = 2 n 2 ). 
The boundary spacetime of a k = 0 topologically-Taub-NUT-AdS solution corresponds to the notrivial (2+1)- 
dimensional part of the Som-Raychaudhuri spacetime m- The value to 2 = 4H 2 (l 2 = 4 n 2 ) corresponds to 
a rotating AdS 3 spacetime, with the four dimensional generalization known as the Rebougas-Tiomno space- 
time jS’O] . 
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The (2 + l)-dimensional metric 1371) contains all the interesting effects exhibited by a four-dimensional 
Godel-type spacetime. We recall that a Godel-type model is the archetypal cosmology exhibiting the prop¬ 
erties associated with the rotation of the universe. Despite the fact that the cosmological rotation is very 
small by the present observation, it cannot be completely ruled out, at least in the early universe m (see 
also m for a up to date discussion and a large set of references). Also, the Godel model is not the first 
but perhaps the best known example of a solution of Einstein’s field equations in which causality may be 
violated and became a paradigm for causality violation in gravitational theory. 

The properties of the Godel-type line elements have been studied by many authors. It was found that, 
for m 2 < 4fl 2 they contain CTC which cannot be removed by going to a covering space and there are no 
reasonable spacelike surfaces. 

Thus asymptotically Taub-NUT-AdS solutions provide an opportunity for the study of quantum field 
theories on Godel-type rotating spacetimes, which is clearly an interesting subject m ■ Given all the problems 
with the Lorentzian interpretation of a Taub-NUT-AdS solution, we should not be surprised if the boundary 
metric and the living CFT present unphysical properties. It is well known M that a generic QFT within 
a causality violating background exhibit a number of pathologies and we may expect to find similar results 
for a (2 + l)-dimensional Godel-type spacetime. 

5 Concluding remarks 

Even if one’s primary interest is in asymptotically flat and AdS solutions, we hope that, by widening the 
context to solutions with NUT charge, one may achieve a deeper appreciation of the theory. There are 
a number of reasons to consider this type of solutions. In some supersymmetric theories, closure under 
duality forces us to consider solutions with NUT parameter. Furthermore, dual mass solutions extremize 
the gravitational action functional in Euclidean quantum gravity and therefore cannot be discarded in spite 
of their causal pathologies. In particular, one may hopes to attain more general features of gravitating 
nonabelian solutions, whether or not containing a magnetic mass. 

Our results seem to invalidate a version of the ”no hair” conjecture for gravitating SU( 2) solutions with a 
magnetic mass, since the structure of these solutions is not completely determined by global charges defined 
at infinity such as mass, NUT parameter, magnetic and electric charges. However, for A = 0, even in the 
presence of a magnetic mass, the decay of the electric potential is too rapid to give a net electric charge 
while w(oo) 2 = 1 corresponds to a zero magnetic charge. These nonabelian solutions are not dyons in the 
usual sense but do posses both electric and magnetic fields. 

As expected, the behavior of the solutions with a negative cosmological constant is different and we find 
configurations with nonzero nonabelian charges. We can attribute this behavior to the different asymptotic 
structure of spacetime. This point was also crucial in obtaining stable hairy black hole solutions in AdS 
spacetime. The presence of a NUT charge does not introduces any new features in this case. Although such 
backgrounds are not of great physical interest, they may provide interesting examples of a more general 
AdS/CFT correspondence. 

We remark also that, on the Lorentzian section, the mass and NUT parameter are unrelated. The 
NUT parameter can be freely specified. An Euclidean section may be obtained by analytically continuing 
the coordinate t and also the parameter n. However, on the Euclidean section, M and n are no longer 
independent quantities EH. Not surprisingly, the nonabelian fields in a Taub-NUT instanton background 
present a very different behavior as compared to the solutions discussed in this paper. 

An important issue that we have not addressed here is the stability of these solutions. Based on the 
experience with n = 0, A < 0 black holes, we expect to find that some of the configurations with no nodes 
in u> are stable solutions against small perturbations. Also, we have all the reasons to expect that all A = 0 
solutions are unstable. The lack of a net electric charge would appear not to improve the chances of stability. 

We can as well to include a positive cosmological constant in the theory and to look for asymptotically 
NUT-de Sitter EYM solutions. Also, a possible extension of this work is to include a dilaton in a theory 
with A = 0. Based on the experience with n = 0 EYM-dilaton solutions EH, we do not expect to find any 
new effects in this case. However, the situation may be different for a supersymmetric theory. An obvious 
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candidate is N = 4, D = 4 gauged SU{2) x SU(2) supergravity, also known as Freedman-Schwarz model 
E53- Nonabelian BPS exact solutions are known in this case for n = 0 and a vanishing electric part of the 
YM potential m It would be interesting to find the corresponding BPS solutions with NUT charge. 

Another interesting issue would be to discuss solutions of a spontaneously broken Einstein-Yang-Mills- 
Higgs theory in the presence of a NUT charge. It can be proven that in this case a magnetic monopole 
necessarily acquires an electric charge. 

We did not addressed here the question of inner structure of these solutions (for r < 77 J. Given the 
special global properties of the vacuum Taub-NUT solution D2f this is clearly an interesting subject, and 
also the study of the possible effects dual mass could have on the singularity structure in a EYM theory. 

One may also speculate about the relevance of these results in string theory and on the possible existence 
of a duality-like transformation relating ordinary hairy black holes to nonabelian solutions possesing a 
magnetic mass. 

A better understanding of these issues may lead to a deeper understanding of the no hair theorems. 
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Figure Captions 


Figure 1: 

The gauge potentials u>(r) and u(r) and the metric functions m(r) and e ~ 2S ^ r ' 1 are plotted as functions of 
radius for typical topological dyon-black hole solutions with A < 0. 

Figure 2: 

The mass-parameter M, the values of the electric and magnetic potentials at infinity uq, ujq and the electric 
and magnetic charges Qe, Qm are plotted as a function of the value of the magnetic potential on the event 
horizon toh, for topological black hole solutions with fixed A, rh, tti. Here and in Figures 4 and 5, Qe and 
Qm are given in units of F/4 tt. 

Figure 3: 

Typical solutions of the EYM theory with A < 0 and a nonvanishing NUT charge. 

Figure 4: 

The mass-parameter M, the values of the electric and magnetic potentials at infinity uq, lu o and the electric 
and magnetic charges Qe, Qm are plotted as a function of the value of the NUT charge for k = 1,0, —1 
solutions with fixed A, rh, Wh, ui- 

Figure 5: 

The mass-parameter M, the values of the electric and magnetic potentials at infinity u 0 , Wo and electric and 
magnetic charges Qe, Qm are plotted as a function of the parameter U\ on the horizon, for k = 1,0,—1 
solutions with fixed A, rh, u>h, n. 

Figure 6: 

Typical one- and two-nodes solutions of the EYM theory with NUT charge and A = 0. The gauge potentials 
u>(r) and u(r) and the metric functions m(r) and e -25 ^ are plotted as functions of radius for fixed values 
of the NUT charge. 
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20 










Figure 2b. 
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Figure 6b. 
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